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HIGHLIGHTS 


►  We  evaluate  a  combination  of  direct-biogas  SOFC  with  MGT  CHP  system. 

►  We  examine  effects  of  reforming  agent,  Uf,  TIT,  compression  ratio  on  system  performance. 

►  The  results  show  the  optimal  operating  parameters  on  SOFC  and  system  performance. 
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The  combination  of  a  direct-biogas  solid  oxide  fuel  cell  (SOFC)  with  a  micro  gas  turbine  (MGT)  system 
offers  great  potential  as  a  green  decentralized  combined  heat  and  power  (CHP)  system.  To  evaluate  the 
potential  use  of  biogas  as  the  main  source  of  energy  for  a  direct-biogas  SOFC-MGT  hybrid  CHP  system, 
a  sensitivity  analysis  was  conducted  under  diverse  operating  conditions  to  investigate  the  influence  of 
key  operating  parameters  of  the  hybrid  CHP  system  with  the  consideration  of  operational  constraints. 
The  key  parameters  in  this  study  were  SOFC  reforming  agent,  SOFC  fuel  utilization  factor  (Of),  turbine 
inlet  temperature  (TIT),  and  compression  ratio.  The  influence  of  variation  in  operating  parameters  on 
plant  performance  was  evaluated  for  the  overall  system  and  SOFC  efficiencies  as  well  as  the  heat-to- 
power  ratio  (TER),  the  power  ratio  of  MGT  to  SOFC  (Pmgt/Fsofc).  and  the  size  of  the  SOFC  stack.  As 
a  reforming  agent  for  direct-biogas  SOFC,  steam  is  more  preferable  than  a  traditional  air— steam  mixture 
in  terms  of  material  limitations  and  SOFC  efficiencies;  however,  an  air— steam  mixture  with  a  small 
amount  of  air  boosts  the  useful  heat  output  and  electricity  generated  by  an  MGT  without  significantly 
affecting  overall  system  efficiency.  The  increase  in  Uf  improves  the  electrical  power  output  produced  by 
the  SOFC  stack,  but  also  requires  more  fuel  to  be  fed  to  the  burner,  resulting  in  an  increase  in  useful  heat 
energy.  Increasing  the  compression  ratio  improves  the  system  electrical  efficiency  but  lowers  useful  heat 
generation;  nevertheless,  increasing  TIT  decreases  the  system  electrical  efficiency  but  improves  the 
efficiency  of  the  CHP  system.  To  achieve  the  optimum  operating  conditions  of  the  hybrid  CHP  system,  the 
operating  parameters  should  be  determined  based  on  the  desired  energy  outcomes. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Biogas  typically  refers  to  a  gas  produced  by  the  biological 
breakdown  of  organic  matter  in  the  absence  of  oxygen.  Organic 
waste  such  as  sludge  from  municipal  wastewater  treatment  plants, 
kitchen  refuse  from  households  and  restaurants,  and  organic  waste 
from  the  food-processing  industry  can  be  converted  into  a  gaseous 
fuel  called  biogas  (a  mixture  of  methane,  carbon  dioxide,  and  other 
minor  gases).  Biogas  is  a  readily  available  but  underexploited 
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energy  source,  because  its  high  levels  of  carbon  dioxide  hinder  its 
use  in  conventional  power-generation  systems,  resulting  in  rela¬ 
tively  low  electrical  conversion  efficiency.  Solid  oxide  fuel  cells 
(SOFCs)  are  a  promising  solution  to  this  problem,  owing  to  their 
tolerance  of  fuel  contaminants  and  flexibility  of  operation  even 
with  diluted  fuel  mixtures. 

Direct  feeding  of  biogas  to  SOFCs  has  been  proven  to  be  feasible 
for  different  SOFC  configurations  and  materials  by  several  experi¬ 
mental  studies  [1-6].  However,  during  the  operation  of  a  direct- 
biogas  SOFC,  carbon  tends  to  form,  which  gradually  deactivates 
the  anode  catalysts  of  the  system.  Shiratori  et  al.  [6]  performed  an 
experimental  study  of  anode-supported  button-cell  SOFCs  fuelled 
by  biogas  with  an  internal  reforming  mode  of  800  °C.  That  study 
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Nomenclature 

Greek  letters 

AC° 

change  of  standard  Gibbs  free  energy,  kj  mol-1 

EcH4 

activation  energy  of  the  methane  reforming  reaction 

AH 

enthalpy  change,  kj  mol-1 

ex 

specific  exergy  flow,  kj  kg-1 

V 

efficiency 

Ex 

exergy  flow  rate,  kW 

* 

rational  efficiency 

F 

faraday  constant,  96,487  C  mol-1 

HRSG 

heat  recovery  steam  generator 

Subscripts 

HX 

heat  exchanger 

0 

properties  of  the  environment 

I< 

equilibrium  constants 

act 

activation 

LHV 

lower  heating  value,  kj  kmol-1 

air 

air,  air  channel 

Mi 

molecular  weight  of  species  i,  kg  mol-1 

ch 

chemical 

m 

mass  flow  rate,  kg  s-1 

CHP 

combined  heat  and  power  overall  system 

ne 

number  of  electrons  participating  in  the 

cone 

concentration 

electrochemical  reaction 

ele 

electrical 

rii 

mole  flow  rate  of  species  i,  kmol  s-1 

FOX 

methane  full  oxidation 

P 

electrical  power,  kW 

fuel 

gas  mixture  at  the  fuel  channel,  fuel  channel 

P 

total  pressure,  kPa 

in 

inlet  value 

Pi 

partial  pressure  of  species  i,  kPa 

inv 

inverter 

Q 

heat  transfer  rate,  kW 

MGT 

micro  gas  turbine 

R 

universal  gas  constant,  8.31434  J  mol-1  K-1 

OC 

open-circuit 

r 

reaction  rate,  kmol  s-1 

ohm 

ohmic 

T 

temperature,  I< 

out 

outlet  value 

TER 

heat-to-power  ratio 

PEN 

positive-electrolyte-negative  structure 

TIT 

turbine  inlet  temperature,  I< 

ph 

physical 

uf 

fuel  utilization  factor 

shift 

shift  reaction 

V 

voltage,  V 

react 

anode  reaction 

mole  fraction  of  species  i 

SOFC 

solid  oxide  fuel  cell  stack 

SR 

steam  reforming 

sys 

overall  system 

revealed  that  the  use  of  air  in  addition  to  actual  biogas  reduces  the 
risk  of  carbon  formation  and  leads  to  more  stable  operation  without 
deteriorating  cell  voltage  due  to  the  lowering  of  anodic  overvoltage. 

In  parallel  with  experimental  studies,  system  models  have 
attracted  tremendous  interest  in  the  past  decade  as  tools  for 
providing  theoretical  guidance  for  SOFC-based  systems.  Yi  et  al.  [7] 
evaluated  integrated  SOFC  reformer  systems  and  found  that  system 
efficiency  drops  insignificantly  by  around  1.1%  when  biogas  is  used 
instead  of  natural  gas.  Piroonlerkgul  et  al.  [8]  investigated  the 
performance  of  biogas-fed  SOFC  systems  utilizing  different 
reforming  agents  (steam,  air,  and  combined  air/steam)  through 
thermodynamic  analysis  to  determine  the  most  suitable  fuel 
processor.  Steam  is  considered  the  most  suitable  reforming  agent, 
as  steam-fed  SOFCs  provide  a  much  higher  power  density  than  air- 
fed  ones,  although  their  electrical  efficiencies  are  slightly  lower. 
While  adding  steam  to  an  air-fed  SOFC  in  the  case  of  a  co-fed  SOFC 
can  improve  power  density,  it  may  lower  electrical  efficiency 
compared  to  an  air-fed  SOFC.  Farhad  et  al.  [9,10]  performed 
a  sensitivity  analysis  of  three  different  biogas-reforming  processors 
in  SOFC  micro-combined  heat  and  power  (CHP)  systems  using 
a  computer  model.  They  considered  anode  exit  gas  recirculation, 
steam  reforming,  and  partial  oxidation.  In  three-dimensional 
computational  fluid  dynamics  (CFD)  simulations  preformed  by 
Vakouftsi  et  al.  [11  ],  mass,  heat,  and  momentum  transfer  equations 
were  combined  with  chemical  and  electrochemical  phenomena 
within  the  inlet  region  of  a  planar  SOFC  unit-cell  configuration 
directly  fed  with  biogas/steam  mixtures.  When  steam  was  added  to 
biogas,  carbon  deposition  was  prevented  but  the  electrical  effi¬ 
ciency  of  the  system  decreased  owing  to  an  increase  in  voltage 
losses. 

Along  with  the  development  of  SOFC  technology,  efforts  have 
been  made  in  the  past  decade  to  integrate  gas  turbines  (GTs)  with 
SOFCs  to  keep  pace  with  rising  demand  for  highly  efficient  energy 


production  while  also  minimizing  environmental  impact.  SOFC-GT 
plants  with  outputs  of  tens  to  hundreds  of  electrical  kilowatts 
(kWe)  have  been  demonstrated  experimentally  [12-17],  making 
conceptual  SOFC-GT  hybrid  systems  feasible.  Because  biogas  is 
increasingly  being  regarded  as  a  potential  renewable  energy  source 
for  distributed  power  generation,  biogas-fed  SOFC-GT  systems 
appear  to  be  one  of  the  most  promising  alternatives  for  distributed 
power  generation.  Moreover,  extending  this  hybrid  system  to  CFIP 
generation  provides  heat  recovery  from  exhaust,  resulting  in  high 
overall  efficiency. 

In  this  study,  a  sensitivity  analysis  was  performed  to  evaluate 
the  influence  of  the  key  operating  parameters  of  a  direct-biogas 
SOFC-micro  gas  turbine  (MGT)  hybrid  CFIP  system  with  an  elec¬ 
trical  power  output  of  200  kWe.  Energy  and  exergy  analyses  were 
used  to  determine  the  causes  of  exergy  losses  and  identify  areas  in 
need  of  improvement  while  adhering  to  material  thermal 
constraints.  Attention  was  paid  to  the  influence  of  air-steam 
mixtures  as  reforming  agents  on  the  direct  internally  reformed 
SOFC  stack  as  well  as  on  the  SOFC-MGT  hybrid  CHP  plant.  The  other 
key  operating  parameters  considered  in  this  study  were  fuel  utili¬ 
zation  factor  (Hf),  turbine  inlet  temperature  (TIT),  and  compression 
ratio.  The  influence  of  variation  in  operating  parameters  on  plant 
performance  was  evaluated  for  the  overall  system  and  SOFC  effi¬ 
ciencies  as  well  as  the  heat-to-power  ratio  (TER)  and  the  power 
ratio  of  MGT  to  SOFC  (Pmgt/Psofc).  Because  of  the  fact  that  the  SOFC 
stack  is  the  most  expensive  part  in  the  initial  investment  cost,  the 
number  of  cells  required  in  the  SOFC  stack  was  also  taken  into 
consideration. 

2.  System  configuration  and  description 

A  schematic  of  the  direct-biogas  SOFC-MGT  hybrid  CHP  system 
used  in  this  study  is  shown  in  Fig.  1.  The  direct  internal-reforming 
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planar  SOFC  is  capable  of  internal  reformation  of  methane  into 
hydrogen  via  dry  reformation,  partial  oxidation,  and  steam  refor¬ 
mation  reactions.  The  other  peripheral  components  include  an 
MGT,  three  gas-to-gas  heat  exchangers  (HX1,  HX2,  and  HX3),  a  heat 
recovery  steam  generator  (HRSG),  a  burner,  a  water  pump,  fuel 
pump,  and  air  pump.  The  operational  parameters  for  the  MGT  and 
other  peripheral  components  are  described  later. 

In  the  integrated  system,  because  the  SOFC  does  not  operate  at 
100%  fuel  utilization,  a  burner  is  needed  to  combust  excess  and 
additional  fuel  to  elevate  the  TIT  to  a  specified  range  for  optimum 
system  performance.  Then,  the  products  of  the  burner  expand  in 
the  turbine  and  the  exhausted  gas  is  further  utilized  by  FIX1  and 
HRSG.  Before  the  flue  gas  is  released  into  the  environment  at 
atmospheric  pressure,  it  is  cooled  to  373  K,  producing  useful  heat. 
The  SOFC  air  and  fuel  streams  are  preheated  by  the  cathode  and 
anode  off-gases  in  HX2  and  HX3  and  are  heated  up  to  1073  K.  An 
inverter  is  also  used  in  the  system  to  convert  the  DC  power  output 
of  the  SOFC  stack  into  AC  power  output.  To  prevent  carbon  depo¬ 
sition  in  the  SOFC,  biogas  is  mixed  with  air  and/or  steam  before  it 
enters  the  cell.  An  HRSG  is  integrated  at  the  MGT  exhaust  outlet, 
supplying  a  predetermined  amount  of  steam  while  a  small  amount 
of  air  from  an  air  pump  is  mixed  with  the  fuel.  For  all  of  the  HXs  and 
the  HRSG,  a  2%  heat  loss  is  assumed.  Possible  variation  in  pressure 
drop  across  each  component  is  assumed  to  be  2%,  except  in  the 
SOFC,  where  it  is  assumed  to  be  3%. 

3.  System  modeling 

The  integrated  system  components  (e.g.,  SOFC,  MGT,  HXs, 
burner,  HRSG,  and  air  and  fuel  pumps)  were  thermodynamically 
modeled  under  steady-state  conditions.  Mass  and  energy  balances 
were  applied  for  each  component  using  lumped  models,  which 
consider  each  component  as  a  control  volume,  with  the  exception 
of  a  one-dimensional  SOFC  model  based  on  previous  work  [18].  The 
following  assumptions  were  made:  (i)  all  system  components 
operate  under  steady-state  conditions;  (ii)  all  gas  mixtures  behave 
as  ideal  gases;  (iii)  the  biogas  is  clean  and  chemically  stable;  (iv) 
heat  loss  from  the  SOFC  is  negligible  and  those  from  HXs  and  the 
HRSG  are  2%  of  the  heat  transferred;  and  (v)  chemical  reactions  are 
in  equilibrium. 

3.1.  Direct-biogas  SOFC  modeling 

The  direct  internal-reforming  SOFC  model  developed  in  a  recent 
work  by  the  authors  [18]  was  modified  to  include  dry-reforming 


and  partial-oxidation  reactions  owing  to  the  presence  of  CO2  and 
02  in  the  anode.  The  model  is  capable  of  capturing  the  distribution 
of  local  temperatures,  species  concentrations,  current  density,  and 
polarization  losses  in  the  streamwise  direction.  The  electro¬ 
chemical  reaction  is  considered  to  be  attributable  to  hydrogen  only; 
the  electrochemical  fuel  value  of  CO  is  readily  exchanged  for 
hydrogen  by  the  rapid  shift  reaction,  assuming  chemical  equilib¬ 
rium  [19]. 

Dry  reforming  is  perhaps  the  most  interesting  option  for  the 
conversion  of  biogas  because  the  major  constituents  of  biogas  are 
carbon  dioxide  and  methane.  However,  the  dry-reforming  reaction 
is  prone  to  carbon  deposition,  which  in  turn  can  cause  rapid 
catalyst  deactivation  [20,21].  Because  steam  and  dry-reforming 
reactions  are  intensively  endothermic,  it  is  necessary  to  supply 
high-temperature  inlet  streams  to  the  SOFC  stack;  however,  this 
reduces  the  overall  efficiency  of  the  system.  This  problem  can  be 
overcome  by  applying  an  exothermic  partial-oxidation  reaction 
that  utilizes  air  as  the  reforming  agent.  However,  the  nature  of  the 
partial  oxidation  of  methane  is  a  matter  of  debate  when  fuel  gas 
and  oxidant  are  mixed.  Hibino  [22,23]  and  Buergler  [24]  have 
suggested  that  in  the  case  of  fuel  in  a  nickel  anode  channel,  partial 
oxidation  of  methane  takes  place  via  a  two-step  mechanism. 
First,  full  oxidation  of  methane  takes  place,  producing  CO2  and 
H20;  then,  synthesis  gas  is  produced  via  steam  reformation  of 
unreacted  methane  while  the  water-gas  shift  reaction  is  at 
equilibrium. 

Mass  balances  are  formulated  for  each  species  on  the  basis  of  the 
relationship  between  the  local  current  and  changes  in  concentra¬ 
tions.  The  following  equations  summarize  the  reactions  considered 
in  the  fuel  cell  stack. 

Full  oxidation  of  CH4  in  the  oxygen-rich  zone  yielding  CO2 
and  H20: 

Full  oxidation  of  methane  :  CH4  +  202->C02  +  2H20  (1) 

Syngas  can  be  generated  by  reforming  reactions  of  unused  methane 
reacting  with  H20  or  C02: 

Steam  reforming  :  CH4  +  H20->C0  +  3H2  (2) 


Dry  reforming  :  CH4  +  C02  ^2CO  +  2H2  (3) 

The  water-gas  shift  reaction  converts  CO  into  H2,  which  is  more 
easily  electrochemically  oxidized  than  CO: 
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Water  -  gas  shift  reaction  :  CO  +  H20<->C02  +  H2  (4) 

Overall  cell  reaction  :  H2  +  ^02  -»H20  (5) 

The  respective  reaction  rate  equation  for  total  oxidation  [25] 
described  in  (Eq.  (1))  on  supported  Ni  catalysts  is  given  as  follows: 


^oc 


A  G° 
~2F~ 


RT 
"2F  ' 


In 


,  ph2°  / 


(12) 


where  the  Faraday  constant  (F)  =  96,485  A  mol-1. 

The  terminal  voltage  of  the  single  cell  can  be  obtained  as 
follows: 


rFOX  = 


fca-PavPo2 


(l  +  'Pch4  +  K0°2X  'P02) 

^ _ ^b'PCH4  Po2 _ 

(i+^c°hVpch4+^2x-Po2 


(6) 


where  pz  is  the  partial  pressure  of  component  i  in  the  gas  mixture. 

The  Arrhenius  reaction  rate  for  component  i  in  combustion 
reaction  lq  is  expressed  as 


Iq  =  A(kk)exp 


RT 


(7) 


V  =  V/0C  -  Vact  -  v0hm  -  Vconc  (13) 

The  terminal  voltage  of  the  cell  is  lower  than  the  open-circuit 
voltage  when  the  current  is  drawn  from  the  cell,  owing  to  three 
types  of  overpotential  losses:  activation  (Vact),  ohmic  (Vohm). 
and  concentration  (VCOnc)  losses.  The  evaluations  of  the  losses 
in  the  present  study  are  the  same  as  those  in  our  previous  report 
[18]  and  are  not  described  here.  The  power  output  of  the  SOFC 
stack  is 

^SOFC  =  Vinv  •  y  [2F(nH2, react  +  nCO, react  +  4nCH4, react)]  H4) 

where  pinv  is  the  inverter  efficiency. 


The  adsorption  constant  for  component  i  in  combustion  reaction 

KPX  is  expressed  as  3.2.  System  parameters  and  efficiencies 


I<PX  =  A(Ki)ex  p(^)  (8) 

The  reaction  rate  of  the  steam-reforming  reaction  is  expressed  as 


rSR 


I<ch4 


Pco-Ph2  \ 

PCH4  "Ph20  '^eq,SR  J 


P  cH4exp 


(9) 


where  FCh4  (=82  kj  mol  ^  is  the  activation  energy  of  the  reaction 
and  /<ch4  (=4274  mol  m~2  bar-1  s-1)  is  the  pre-exponential  factor 
(Table  1). 

The  dry-reforming  reaction  (Eq.  (3))  is  not  explicitly  included 
but  is  implicitly  considered  through  the  water-gas  shift  reaction 
(Eq.  (4)).  As  mentioned  above,  the  water-gas  shift  reaction  was 
always  assumed  to  be  at  equilibrium  in  this  study.  The  equilibrium 
constant  for  this  reaction  and  the  steam-reforming  reaction  can  be 
described  as  follows  [26]: 


Kp shift  =  exp^  -  0.2935Z3  +  0.6351Z2  +  4.1788Z  +  0.3619^ 

(10) 


Keq.SR  =  1.0267  x  1010exp(  -  0.2513Z4  +  0.3665Z3 

+  0.5810Z2  -27.134Z  +  3.2770)  (11) 

where  Z  =  (1000/T)  -  1. 

The  open-circuit  voltage  (Voc)  is  described  by  the  Nernst 
equation  as  a  function  of  operating  temperature  (T)  and  partial 
pressure  (p): 


Because  the  fuel  is  not  completely  consumed  by  electro¬ 
chemical  reactions  in  the  anode  channel,  the  rate  of  excess  fuel 
flow  is  determined  by  Uf,  a  key  parameter  investigated  in  this 
study.  Uf  is  defined  here  as  the  ratio  of  fuel  consumed  by 
anode  reactions  to  the  fuel  entering  anode  channels,  and  is 
expressed  as 


Uf 


nH2,  react  +  ^CO, react  +  4^CH4, react 
(nH2,SOFC  +  ftCO,SOFC  +  4nCH4,SOFc)in 


(15) 


The  TER  of  the  combined  plant  indicates  the  ratio  of  useful 
thermal  energy  to  electricity  generation.  It  is  determined  as 


TER 


Quseful 


Fsofc  +  Fmgt 


(16) 


We  defined  two  types  of  efficiencies:  energetic  (or  fuel)  effi¬ 
ciency  and  rational  efficiency  for  steady-state  processes,  as 
described  by  Kotas  [27]. 


3.2.1.  Energetic  efficiency  (rj) 

In  any  system,  energetic  efficiency  is  defined  as  the  ratio  of 
energy  in  the  product  output  to  the  energy  in  the  fuel  input.  It  can 
be  applied  practically  as  follows: 

The  SOFC  efficiency  can  be  determined  as 


^SOFC  — 


_ Fsofc _ 

S  (^anode  ’  LH V)  jn 

ch4,co,h2 


X  100% 


(17) 


The  system  electrical  efficiency  for  a  power  generation  system  is 
defined  as 


Table  1 

Calculation-based  parameters  for  the  reaction  rate  constants  in  Eq.  (7)  and  the  adsorption  constants  in  Eq.  (8)  (taken  from  Ref.  [25]). 


A(kk) 

Ek  (1<J  mol  1) 

A(IQ 

AHfe  (1<J  mol  b 

/<a(l<mol  kgcajbar  2h  !) 
fcb(kmol  kg“tbar 2h 1) 

2.92  x  106 

86.0 

^(bar1) 

1.26  x  10"1 

-27.3 

2.46  x  106 

86.0 

/COXfbai-1) 

7.87  x  10~7 

-92.8 
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Ve\e,sys  ~ 


[  SOFC  ■ 


!  MGT 


gas '  LHVbiogas 


x  100% 


(mbiog< 

The  CHP  system  efficiency  is  defined  as 


^CHP 


'  SOFC  ■ 


1  MGT  ■ 


'  Quseful 


(^biog, 


as  LHVbiogas 


x  100% 


(18) 


(19) 


To  gain  more  insight  into  the  system,  exergy  analysis  was  used 
to  identify  and  locate  irreversibility.  The  exergy  associated  with 
a  material  stream  is  equal  to  the  maximum  amount  of  work 
obtainable  when  the  stream  is  brought  from  its  initial  state  to  the 
dead  state.  The  exergy  transfer  rate  associated  with  a  material 
stream  can  be  divided  into  physical  and  chemical  exergy 
components: 

Ex  =  n(exph  +  exch)  (20) 

Physical  exergy  is  the  work  obtainable  by  taking  the  substance 
through  reversible  processes  from  its  initial  state  at  temperature  T 
and  pressure  P  to  the  environmental  state.  It  can  be  calculated  as 

exph  =  -  ho)  ~  T0(s  -  So)],-  (21) 

i 

where  h  is  the  specific  enthalpy,  s  is  the  specific  entropy,  and  the 
properties  indicated  with  the  subscript  0  refer  to  the  environ¬ 
mental  state. 

Chemical  exergy  is  equal  to  the  maximum  amount  of  work 
obtainable  when  the  substance  under  consideration  is  brought 
from  the  environmental  state,  defined  by  the  parameters  To  and  Po, 
to  the  reference  state  by  processes  involving  heat  transfer  and 
exchange  of  substances  only  with  the  environment.  The  chemical 
exergy  for  mixtures  can  be  expressed  as 


exch  =  Ylx‘'exo>  +  RT°Ylx‘Mx‘)  (22) 

i  i 

where  ex0,j  denotes  the  standard  molar  chemical  exergy  of  the  i-th 
species,  assuming  the  reference  atmospheric  composition  given  by 
Kotas  [27].  The  exergy  transfer  connected  with  the  heat  transfer 
rate  Q.  at  temperature  T  can  be  expressed  as 


EXq 


(23) 


3.2.2.  Rational  efficiency  (\p) 

Rational  efficiency  is  defined  by  Kotas  [27]  as  the  ratio  of  the 
desired  exergy  output  to  the  exergy  used  or  consumed.  It  provides 
a  realistic  indicator  of  a  system’s  efficiency  by  applying  the  exergy 
of  the  working  fluid  at  the  outlet  of  a  component  to  the  next 
component. 

The  SOFC  rational  efficiency  is  defined  as 


fc)FC  — 


^SOFC 


(Ex, 


fuel  ' 


Ex, 


SOFC, in 


(EXfuei 


+  Ex 


X  100% 


SOFC, out 


(24) 


The  rational  efficiency  for  CHP  system  is  defined  as 


,  ^SOFC  +  +  EXq 

YCHP  =  - —  X  100% 


EXfuei  ;SyS  Ex( 


exhaust 


(25) 


Table  2 

Operational  parameter  values  for  the  SOFC  simulation. 


Parameters 

Value 

Biogas 

CH4  60%:C02  40% 

Stack  input  data 

Fuel  utilization  factor  Uf(-) 

0.75 

Air  utilization  factor  Ua  (-) 

0.25 

Average  current  density  (A/m2) 

4000 

Number  of  channels  per  SOFC  cell 

20 

Air  inlet  temperature  to  the  SOFC  Tair,in  (K) 

1073 

Fuel  inlet  temperature  to  the  SOFC  TfUei,in  (K) 

1073 

Cell  length  (mm) 

100.0 

Width  covered  by  one  channel  (mm) 

5.0 

Air  channel  height  (mm) 

1.5 

Fuel  channel  height  (mm) 

0.4 

Anode  thickness  (pm) 

500 

Cathode  thickness  (mm) 

50 

Electrolyte  thickness  (pm) 

10 

4.  Simulation  results 

4.1.  Computational  conditions 

Because  the  SOFC  stack  is  considered  the  central  part  of  the 
hybrid  system,  the  simulation  results  are  presented  in  two  parts. 
First,  the  performance  of  the  biogas-fed  SOFC  operated  at  atmo¬ 
spheric  pressure  utilizing  different  reforming  agents  (steam  and 
combined  air/steam)  was  investigated  via  thermodynamic  analysis 

to  determine  the  most  suitable  feed,  using 

the  operational 

parameter  values  listed  in  Table  2.  SOFC  performance  under  co-flow 
operation  was  also  analyzed  because  it  generally  has  a  more 
uniform  temperature  distribution  than  other  flow  configurations 
[18].  These  results  are  discussed  first.  The  second  part  presents  the 
sensitivity  analysis  of  the  direct-biogas  SOFC-MGT  hybrid  CHP 
system  shown  schematically  in  Fig.  1.  To  understand  the  opera- 

tional  scenarios  of  the  hybrid  CHP  system, 

an  independent 

parameter  analysis  of  a  single  component  is  not  enough  to  assess 
the  whole  system,  because  all  components  in  the  system  affect  one 
another.  The  input  operational  parameter  values  for  the  SOFC  and 
the  other  system  components  presented  in  Tables  2  and  3  were 

used  as  constants  throughout  the  study, 
otherwise. 

unless  indicated 

The  inlet  gas  compositions  are  determined  by  considering  the 
carbon-hydrogen-oxygen  (C— H— O)  ternary  diagrams  shown  in 
Fig.  2.  In  principle,  it  is  feasible  to  directly  feed  biogas  containing 
the  natural  reforming  agent,  CO2,  into  a  high-temperature  SOFC 
without  an  additional  reforming  agent.  However,  in  regard  to  the 
most  common  form  of  biogas  considered  in  this  work  (represented 

Table  3 

Operational  parameter  values  for  system  simulation. 

Parameters 

Value 

Total  electrical  net  output  power  (kW) 

200 

Steam-to-carbon  ratio  (-) 

2 

Compression  ratio  (-) 

8 

Turbine  inlet  temperature  TIT  (I<) 

1073 

Exhaust  gas  temperature  (I<) 

373 

Pump  isentropic  efficiency  7]pump  (-) 

0.95 

Air  compressor  isentropic  efficiency  (-) 

0.75 

Fuel  compressor  isentropic  efficiency  (-) 

0.75 

Turbine  isentropic  efficiency  (-) 

0.85 

Pump  mechanical  efficiency  (-) 

0.98 

Air  compressor  mechanical  efficiency  (-) 

0.98 

Fuel  compressor  mechanical  efficiency  (-) 

0.98 

Turbine  mechanical  efficiency  (-) 

0.98 

Inverter  efficiency  ^inv  (-) 

0.95 
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Fig.  2.  C— H— 0  ternary  diagram  with  a  carbon  deposition  boundary  at  1073,  973,  and 
873  K  and  1  atm. 

by  60%  CH4  and  40%  CO2  in  volume),  biogas  lies  above  the  carbon 
deposition  boundary  curves,  indicating  that  solid  carbon  exists  in 
heterogeneous  equilibrium.  The  location  of  biogas  may  be  moved 
below  the  carbon  deposition  boundary  by  adding  steam  or  oxygen. 
As  clearly  shown  in  Fig.  2,  an  increase  in  steam  or  oxygen  as 
a  reforming  agent  for  methane  in  biogas  can  minimize  the  risk  of 
carbon  deposition.  In  this  study,  six  different  inlet  gas  compositions 
were  examined.  Cases  I  and  VI  are  shown  in  Fig.  2;  the  others  (cases 
II— V)  are  not  shown  in  the  figure  but  are  located  somewhere 
between  these  two  points.  Detailed  gas  compositions  are  explained 
in  the  next  section  (Fig.  3). 

4.2.  Direct-biogas  SOFC  simulation 

A  strong  endothermic  reaction  due  to  the  presence  of  steam  can 
lead  to  local  temperature  gradients,  especially  near  the  entrance  of 
the  stack,  resulting  in  mechanical  failure  due  to  thermally  induced 
stress.  In  this  study,  the  maximum  allowable  temperature  gradient 
and  the  maximum  allowable  cell  temperature,  which  are  the  most 
important  operational  constraints  for  a  planar  SOFC,  were  set  to 
1300  I<  and  5  I<  mm-1,  respectively,  following  Stiller  et  al.  [28]. 


Fig.  3.  SOFC  cell  temperature  profiles  for  different  air-steam  mixtures  as  SOFC 
reforming  agents. 


Table  4  presents  the  sensitivity  of  the  system  to  different  ratios 
of  air  and  steam  to  fuel  for  cases  I— VI,  illustrating  how  reforming 
agents  may  affect  SOFC  performance.  As  mentioned  in  the  previous 
section,  anode  feed  gases  were  located  below  the  carbon  deposition 
boundary  in  all  cases.  As  shown  in  Table  4,  an  increase  in  air  input 
into  the  anode  deteriorated  the  SOFC  energetic  and  rational  effi¬ 
ciencies  from  50.8%  to  47.4%  and  from  72.3%  to  70.0%,  respectively, 
mainly  due  to  partial  oxidation.  The  SOFC  cell  temperature  profiles 
for  the  six  different  ratios  of  air  and  steam  to  fuel  are  depicted  in 
Fig.  2.  The  presence  of  oxygen  caused  the  cell  temperature  near  the 
channel  inlets  to  rise,  and  consequently  accelerated  the  strong 
endothermic  reactions  that  take  place  during  the  steam  reforma¬ 
tion  of  methane.  As  the  steam-reforming  reaction  rate  increases, 
the  cell  temperature  drops  more  rapidly  and  the  cell  temperature 
gradient  becomes  steeper  near  the  gas  inlets.  The  maximum  cell 
temperatures  (TPEn,  max)  and  the  maximum  cell  temperature 
gradients  (5TpEN/9*)max  corresponding  to  Fig.  2  are  listed  in  Table  4. 
According  to  the  table,  only  cases  I— IV  led  to  safe  operation  of  the 
system  under  the  material  constraints.  The  most  favorable  oper¬ 
ating  conditions  were  in  case  I,  which  had  the  lowest  maximum 
temperature  gradient  and  the  highest  energetic  (ijsofc)  and  rational 
efficiencies  (^ sofc)-  It  should  be  noted,  however,  that  a  high  exhaust 
gas  temperature  is  favorable  for  hybrid  operation.  The  effects  of 
reforming  agents  on  the  hybrid  CFIP  system  are  discussed  later  in 
this  report. 

4.3.  System  simulation 

A  sensitivity  analysis  was  used  to  quantify  the  effects  of  air- 
steam  mixtures  as  reforming  agents,  Df,  TIT,  and  the  compression 
ratio  on  system  performance  and  the  size  of  the  SOFC  stack. 

4.3 A.  Influence  of  air-steam  mixtures  as  reforming  agents 

The  efficiencies  of  the  direct-biogas  SOFC-MGT  hybrid  CHP 
system  under  various  air-steam  mixtures  as  reforming  agents 
(anode  feed  gas  compositions  of  cases  I-V  in  Table  3)  are  depicted 
in  Fig.  4.  As  the  amount  of  air  input  increased,  both  the  SOFC 
energetic  (tjsofc)  and  the  SOFC  rational  efficiencies  (i/'sofc) 
decreased,  while  the  exhaust  heat  from  the  SOFC  increased.  This 
compensated  for  the  drop  in  electrical  power  produced  by  the  SOFC 
stack.  These  results  suggest  that  air-steam  mixtures  have  only 
slight  effects  on  system  electrical  efficiency  (^eie,syn)»  CFIP  system 
efficiency  (??chp),  and  CFIP  system  rational  efficiency  (^chp)-  The 
effects  of  air-steam  mixtures  as  reforming  agents  on  the  required 
number  of  SOFC  cells,  the  output  power  ratio  of  MGT  to  SOFC  (Pmgt/ 
Psofc),  and  the  TER  of  the  direct-biogas  SOFC  in  the  hybrid  CHP 
system  are  shown  in  Fig.  5.  At  an  air:biogas  ratio  of  0.4,  the  required 
number  of  cells  decreased  by  approximately  4%,  whereas  Pmgt/ 
Psofc  increased  from  0.37  to  0.40.  As  more  air  was  added,  TER 
increased  from  0.38  to  0.41  owing  to  more  heat  energy  in  the  SOFC 
off-gases  going  to  the  MGT  system,  leading  to  a  reduction  in  the  fuel 
fed  to  the  burner.  It  should  be  noted  that  only  the  anode  feed  gas 
composition  of  case  I  was  used  for  the  rest  of  the  system  evaluation. 


Table  4 

Summary  of  SOFC  performance  at  atmospheric  pressure  based  on  different  SOFC 
reforming  agents. 


Case 

Biogas:  steam:  air 

VSOFC  (%) 

u 

0 

fpEN,  max  (K) 

(dTpEN/dX  )  max(  k/  111  m  ) 

I 

1:2.0:0.0 

50.8 

72.3 

1077.1 

1.6 

II 

1:1.9:0.1 

50.2 

71.9 

1081.5 

2.4 

III 

1:1.8:0.2 

49.6 

71.5 

1086.3 

3.3 

IV 

1:1.7:0.3 

48.9 

71.2 

1091.1 

4.3 

V 

1:1.6:0.4 

48.0 

70.7 

1096.2 

5.4 

VI 

1:1.5:0.5 

47.4 

70.0 

1101.7 

6.5 
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Fig.  6.  Influence  of  the  fuel  utilization  factor  on  SOFC  and  system  efficiencies. 


Fig.  4.  Influence  of  air-steam  mixtures  as  reforming  agents  on  SOFC  and  system 
efficiencies. 


4.3.2.  Influence  of  Uf 

The  Uf  of  a  commercial  SOFC  is  75-85%  [29].  Uf  is  an  important 
parameter,  which  is  closely  related  to  the  performance  of  SOFCs 
and  exhaust  heat,  reflecting  the  performance  of  hybrid  CFIP 
systems.  The  proper  amounts  of  exhaust  heat  and  unreacted  fuel 
from  the  electricity  conversion  process  result  in  an  efficient  GT 
cycle.  The  influence  of  Uf  on  SOFC  and  system  performance  is 
shown  in  Fig.  6.  As  Uf  increased  from  0.65  to  0.90,  SOFC  energetic 
efficiency  (i^sofc)  and  CFIP  system  efficiency  (tjchp)  increased  from 
42.3%  to  59.8%  and  from  67.1  to  69.7%,  respectively.  The  significant 
increase  of  psofc  is  due  to  the  fact  that  as  Uf  is  enhanced,  electro¬ 
chemical  reaction  rates  increase,  thus  raising  the  electrical  power 
output  produced  by  the  SOFC  stack.  However,  when  considering 
exergy  in  the  outgoing  streams,  the  SOFC  rational  efficiency  (i/'sofc) 
peaked  at  76.2%  when  Uf  was  0.65,  while  CHP  system  rational 
efficiency  (i/'chp)  was  relatively  constant  around  62%.  The  influence 
of  Uf  on  the  required  number  of  SOFC  cells,  Pmgt/Fsofg  and  TER  is 
presented  in  Fig.  7.  At  a  Uf  of  0.75,  the  required  number  of  cells  was 
lowest  (5277)  among  all  of  the  cases  studied.  This  is  mainly  due  to 
compensation  between  the  electrical  power  produced  by  the  MGT 


and  the  SOFC.  As  Uf  changed  from  0.65  to  0.90,  Pmgt/Fsofc  dropped 
from  0.41  to  0.36,  owing  to  a  considerable  increase  in  the  SOFC 
efficiency.  When  Uf  increases,  it  results  in  less  unused  fuel  in  the 
anode  exhaust,  and  consequently,  more  additional  fuel  is  needed  to 
maintain  the  set  TIT.  At  the  same  time,  TER  increased  from  0.34  to 
0.40,  because  the  amount  of  heat  in  flue  gas  increased  as  the 
amount  of  additional  fuel  fed  to  the  burner  increased. 

4.3.3.  Influence  of  compression  ratio  and  TIT 

The  compression  ratio  and  TIT  are  key  design  parameters  of  an 
MGT-SOFC  hybrid  system,  because  they  considerably  affect  the  heat 
balance  between  the  MGT  and  SOFC  units.  In  the  present  study,  the 
compression  ratio  and  TIT  were  varied  from  2  to  12  and  973  to 
1273  K,  respectively.  In  Fig.  8,  the  effects  of  compression  ratio  and 
TIT  on  the  CHP  system  efficiency  (^chp)  and  CHP  system  rational 
efficiency  (i/'chp)  are  plotted.  It  should  be  noted  that  when  the 
compression  ratio  is  set  at  12  with  an  input  operational  Uf  of  0.75, 
TIT  exceeds  973  K.  Therefore,  no  data  point  for  this  case  is  shown  in 
the  figure.  As  can  be  seen  in  the  figure,  at  a  compression  ratio 
between  2  and  3,  the  deviation  of  tjc hp  was  not  significant. 
However,  when  the  ratio  was  increased  from  3  to  12,  tjchp 
decreased  linearly.  This  is  because  increasing  the  compression  ratio 
reduces  the  amount  of  useful  heat  available,  resulting  in  a  decrease 
in  ?7chp-  On  the  other  hand,  when  considering  exergy  in  exhaust  gas, 
i/'chp  increased  with  an  increase  in  the  compression  ratio,  and  the 
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Fig.  5.  Influence  of  air-steam  mixtures  as  reforming  agents  on  Pmgt/Psofc.  TER,  and  Fig.  7.  Influence  of  the  fuel  utilization  factor  on  Pmgt/Psofc.  TER,  and  the  required 
the  required  number  of  SOFC  cells.  number  of  SOFC  cells. 
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□  TIT  973  K  o  TIT  1073  K  Aim  173  K  v  TIT  1273  K 


□  TIT  973  K  o  TIT  1073  K  a  TIT  1173  K  v  TIT  1273  K 


optimum  i/'chp  was  achieved  at  a  compression  ratio  in  the  range  of 
9-11.  In  addition,  increasing  TIT  considerably  improved  the  ijchp 
and  ^chp  of  the  system.  This  is  because  the  production  of 
byproduct  heat  is  significantly  boosted  with  increasing  TIT. 

□  TIT  973  K  o  TIT  1073  K  a  TIT  1173  K  v  TIT  1273  K 


□  TIT  973  K  o  TIT  1073  K  a  TIT  1 173  K  v  TIT  1273  K 


However,  as  TIT  increases,  more  additional  fuel  is  supplied  to  the 
burner,  leading  to  an  increase  in  the  system  fuel  consumption 
rate.  This  also  results  in  lower  net  system  electrical  efficiency. 
Indeed,  as  can  be  seen  in  Fig.  9,  increasing  TIT  in  the  present 
study  had  a  negative  impact  on  system  electrical  efficiency  (rye le.sys) 
due  to  an  increase  in  fuel  feed.  In  addition,  increasing  the 
compression  ratio  from  2  to  12  enhanced  ^eie.sys  by  approximately 
18%  (Fig.  9).  The  system  electrical  efficiency  increased  more  rapidly 
at  a  low  compression  ratio.  As  shown  in  Fig.  10,  the  TER  increased 
with  an  increase  in  TIT,  as  more  heat  that  is  useful  was  produced. 
Nevertheless,  the  TER  decreased  with  an  increase  in  compression 
ratio,  because  the  MGT  gains  efficiency  as  electricity  is  generated. 
Variation  in  Pmgt/Psofc  according  to  compression  ratio  and  TIT  is 
shown  in  Fig.  11.  The  Pmgt/Psofc  reached  a  maximum  value  at 
a  compression  ratio  between  5  and  10  when  TIT  ranged  from  973  to 
1273  K.  The  maximum  value  of  Pmgt/Psofc  shifts  toward  a  higher 
compression  ratio  as  TIT  increases.  As  can  be  seen  in  Fig.  12,  at 
lower  compression  ratios,  the  required  number  of  SOFC  cells 
decreased  significantly  and  reached  a  minimum  value  at 
a  compression  ratio  in  the  vicinity  of  8.  In  addition,  increasing  TIT 
reduced  the  required  number  of  SOFC  cells.  This  is  mainly  because 
the  portion  of  electrical  power  generated  by  MGT  increases  as  TIT 
increases,  leading  to  a  reduction  in  the  required  number  of  SOFC 
cells. 


□  TIT  973  K  o  TIT  1073  K  a  TIT  1173  K  v  TIT  1273  K 


Fig.  12.  Influence  of  compression  ratio  and  TIT  on  the  required  number  of  SOFC  cells. 
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5.  Conclusions 

Performance  evaluations  of  a  direct-biogas  SOFC-MGT  hybrid 
CHP  system  were  conducted.  It  is  important  to  closely  monitor  the 
operational  performance  of  SOFCs  to  minimize  SOFC  degradation 
due  to  thermal  constraints  and  carbon  deposition  on  Ni-based 
anodes.  The  main  SOFC  model  used  in  this  study  was  based  on 
previous  work  [18],  whereas  other  system  components  were  zero¬ 
dimensional  component  models.  Sensitivity  analysis  was  carried 
out  to  investigate  the  influences  of  the  main  variables  on  the 
system.  The  main  parameters  considered  were  air— steam  mixtures 
as  reforming  agents,  Uf,  TIT,  and  the  compression  ratio.  Based  on  the 
results  of  our  simulations,  the  following  conclusions  can  be  made: 

•  Considering  individual  direct-biogas  SOFC  operation,  as  the 
average  cell  operating  temperature  increases  as  more  air  is 
added  to  the  biogas,  carbon  deposition  is  less  likely  to  form  in 
the  anode  channels.  Nevertheless,  the  presence  of  oxygen  has 
a  negative  impact  on  SOFC  performance  and  also  causes 
temperature  stress  near  the  stack  inlet  due  to  partial  oxidation. 

•  The  addition  of  a  small  amount  of  air  to  biogas  does  not  have 
a  significant  effect  on  ^eie,syn*  ?7chb  or  i/'chp  of  the  direct-biogas 
SOFC  in  the  hybrid  CFIP  system.  Flowever,  it  raises  the  stack 
operating  temperature  and  leads  to  an  increase  in  useful  heat 
output  as  well  as  electrical  power  production  by  the  MGT, 
which  in  turn  reduces  the  number  of  SOFC  cells  required. 

•  In  the  system  studied,  the  electrical  power  output  produced  by 
the  SOFC  stack  was  directly  proportional  to  Uf,  however,  ^eie.syn 
was  not  significantly  affected  by  variation  in  Uf.  The  smallest 
number  of  cells  was  achieved  at  a  Uf  of  0.75,  whereas  the  TER 
increases  with  an  increase  in  Uf. 

•  Increasing  the  compression  ratio  improves  ^eie.sys  but  reduces 
TER,  whereas  increasing  TIT  has  the  opposite  outcomes.  Flowever, 
increasing  the  compression  ratio  and  TIT  has  the  same  influence 
on  \f/Cm  Fmgt/Psofg  and  the  required  number  of  SOFC  cells. 
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